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^ ■ Abstract 

The multicomponent 2D Toda hierarchy is analyzed through a factorization problem associated to an infinite- 
dimensional group. A new set of discrete flows is considered and the corresponding Lax and Zakharov-Shabat 
Cn| ■ equations are characterized. Reductions of block Toeplitz and Hankel bi-infinite matrix types are proposed and 
studied. Orlov-Schulman operators, string equations and additional symmetries (discrete and continuous) are 
considered. The continuous-discrete Lax equations are shown to be equivalent to a factorization problem as well 
' as to a set of string equations. A congruence method to derive site independent equations is presented and used 
■ to derive equations in the discrete multicomponent KP sector (and also for its modification) of the theory as 
well as dispersive Whitham equations. 

^ 1 Introduction 

> . . . . 

' This paper revisits the multicomponent 2D Toda hierarchy [30j from the point of view of the factorization problem 
^SJ , associated to an infinite-dimensional group. Our main motivation is the recent discovery [3] of underlying 
' integrable structures of multicomponent type in the theory of multiple orthogonal polynomials which is in turn 
, connected to models of non-intersecting Brownian motions. Having in mind the fruitful applications of the Toda 
0^ ' hierarchy to the theory of orthogonal polynomials and to the Hermitian random matrix model (see for instance 
. |14|-[21j). it is expected that the formalism of multicomponent integrable hierarchies can be similarly applied to 
the study and characterization of multiple orthogonal polynomials and non-intersecting Brownian motions. In 
particular, the semiclassical (dispersionless) limit of multicomponent integrable hierarchies should be relevant 
for the analysis of large N) type limits, see for instance [22j. An important piece of the technique required for 
these applications was recently provided by Takasaki and Takebe [571 US]. Indeed, they proved that the universal 
^ • Whitham hierarchy (genus case) [161 can be obtained as a particular dispersionless limit of the multicomponent 
5^ I KP hierarchy. 

The applications of the Toda hierarchy to the characterization of semiclassical limits make an essential use of 
the notion of string equations [14j-|21|-[10|. In recent years the formalism of string equations for dispersionless 
integrable hierarchies '26] has been much developed [321 119 1 but, to our knowledge, a similar formalism for 
dispersive multicomponent integrable hierarchies is not yet available. One of the main goals of this paper is to 
extend the formalism of string equations to multicomponent 2D Toda hierarchies. In this sense the consideration 
of factorization problems for these hierarchies turns to be of great help in order to introduce basic ingredients 
such as discrete flows, Orlov-Schulman operators and additional symmetries. 

The theory of the multicomponent KP hierarchy is discussed in length in the papers [Hill], see also [20] for its 
applications to geometric nets of conjugate type. In [30j it was noticed that r functions of a 2A'^-multicomponent 
KP provide solutions of the iV-multicomponent Toda hierarchy. The introduction of integer parameters in 
the multicomponent KP hierarchy goes back to [5] and the corresponding discrete flows, which are used in 
two different ways in [TS] [4] , are essential for the derivation of the dispersionless Whitham hierarchy from the 
multicomponent KP hierarchy [571 HE] , [IS] ■ In the present paper we introduce a set of discrete flows for the 
multicomponent 2D Toda hierarchy. Its role in the formulation of the corresponding dispersionless limits will 
be discussed in length in a forthcoming paper. 
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The layout of the paper is as follows: In §1 we introduce a factorization problem in a Lie group as the one 
presented in [5D] . This factorization problem is rooted in the ideas used for the KP case in [531 U j for the so 
called discrete KP hierarchy. Then we derive the continuous and discrete Lax equations for the multicomponcnt 
2D Toda hierarchy. We notice that in our discussion the set of discrete flows, which to our knowledge where 
not considered before for this hierarchy, are formulated in equal footing to the continuous flows. We also 
show some examples of members of the hierarchy and, in particular, multicomponent equations of Toda type 
involving partial difference operators only, or combined partial difference and partial derivatives. We end this 
section with the formulation of several classes of reductions of the multicomponent 2D Toda hierarchy involving 
biinfinite block Toeplitz and Hankel matrices. The consideration of these types of reductions is motivated by 
their relevance in integrable hierarchies such as the infinite Toda hierarchy [2j or the Ablowitz-Ladik lattice 
hierarchy [6] . For some reductions we characterize solutions of the hierarchy which are periodic in the discrete 
variables. Moreover, for the Hankel case we get generalizations of the bigraded reduction, see [7], associated 
with extended flows of the 1-component ID Toda hierarchy [fi . 

In §2 we formulate the theory of string equations for the multicomponcnt 2D Toda hierarchy. We start by 
defining the Orlov-Schulman operator ^23 and then we derive its Lax equations from the factorization problem 
introduced in §1. We also show how the Lax equations imply in turn the factorization problem. In this way 
we stablish the equivalence between the factorization problem and the extended Lax formulation, involving 
discrete flows and the Orlov-Schulman operator, of the multicomponent 2D Toda hierarchy. Moreover, we also 
prove the equivalence between the extended Lax formulation and a particular type of string equations for the 
multicomponent 2D Toda hierarchy. This generalizes the result for the one-component case stablished in [26] . 
Finally, we use the factorization problem and the canonical pair of Lax and Orlov-Schulman operators to provide 
a natural formulation of the additional symmetries of the multicomponent 2D Toda hierarchy. As a consequence 
we characterize the string equations as invariance conditions under additional symmetries. 

The paper ends with two appendices. In the first appendix the congruence method for deriving n-independent 
equations is shown. It is applied to get the main equations of the discrete multicomponcnt KP hierarchy: N-wave 
equations, Darboux equations and multiquadrilateral lattice equations [9j . We also use this method to formulate 
the dispersive Whitham equations in terms of scalar Lax and Orlov-Schulman opeartors, which constitute the 
starting point for the discussion of the dispersionless limits of the multicomponcnt 2D Toda hierarchy. Finally, 
the second appendix contains the proofs of the main Propositions of the paper. 

1.1 Lie algebra setting 

In this paper we only consider formal series expansions in the Lie group theoretic setup without any assumption 
on their convergency. We also remark that along the paper we use the following notations. For given Lie algebras 
01 C 02, and AT, y e 02 then X ~ Y + gi means X — Y E gi. For any Lie groups Gi C G2 and a,b E G2 then 
a = G2 ■ b stands for a ■ G G2. Let {Eki}^i^i be the standard basis {Eki)k'V = ^/i' ^fefe' of Mm[C) and Dat 
denote the identity matrix in Mn{£,). We also denote the algebra unit as 1. 

If M7v(C) denotes the associative algebra of complex N x N matrices we will consider the linear space of 
sequences 



The shift operator A acts on these sequences as (A/)(n) :— f{n + 1). A sequence X : Z — > AIn{C) acts by left 
multiplication in this space of sequences, and therefore we may consider expressions of the type XA^ , where 
X = X{n) is a sequence which acts by left multiplication: {XA^){f){n) := X{n) ■ f{n + j). 

Moreover, defining the product (X(n)A') • (F(n)A^) := X{n)Y{n + i)^^^^ and extending it linearly we have 
that the set g of Laurent series in A is an associative algebra, which under the standard commutator is a Lie 
algebra. Observe that g can be thought either as Mz(MAr(C)), i.e. bi- infinite matrices with Mjv(C) coefficients, 
or as Mm{Mj_{£,)), i.e. N x N matrices with coefficients bi-infinite matrices. 

This Lie algebra has the following important splitting 
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are Lie subalgebras of g with trivial intersection. 



1.2 The Lie group and the factorization problem 

The group of linear invertible elements in q will be denoted by G and has q as its Lie algebra, then the splitting 
H]) leads us to consider the following factorization oi g E G 

9 = 9-^-9+, 9±&G± (2) 

where G± have g± as their Lie algebras. Explicitly, G+ is the set of invertible linear operators of the form 
^jyf)gj{n)hP \ while G_ is the set of invertible linear operators of the form 1 + '^j^Q9j{n)KK 
An alternative factorization is the Gauss factorization 

9^9-^-9+, 9±&G± (3) 

where G+ is the set of invertible linear operators of the form go,+ ('^) + X]j>o do i^)^'' ! with 50,+ ■ Z GL(iV, C) 
an invertible upper triangular matrix, while G- is the set of invertible linear operators of the form go^-(ri)~^ + 
Sj<o ilji'^)^^ with 5o,- : 2 GL(7V, C), such that 50,- = Dtv + A, being A a strictly lower triangular matrix in 
Mm{C). If the factorization ([3]) exists then it will also exist ^ by defining g+ = go,- '9+, 5- = 9q,- ■ 5-- The 
elements g with a factorization ([3]) are said to belong to the big cell [4j , hence the factorization can be considered 
only locally. Thus, we will consider elements g in the big cell so that the factorization ([2|) holds, avoiding the 
generation of additional problems connected with these local aspects. 

Now we introduce two sets of indexes, S = {!,..., N} and S = {!,..., N}, of the same cardinality N. In 
what follows we will use letters k, I and k, I to denote elements in S and S, respectively. Furthermore, we will 
use letters a, b, c to denote elements in 5 := S U S. 

We define the following operators Wq, Wq £ G 

N 

Wo:=Y.EkkA''e^T^ot.>^^\ (4) 
fc=i 

N 

Wo := J2 EkkA-'^' e^r=i 'r^^'' (5) 

k=l 

where Sa G Z, tja G C are deformation parameters, that in the sequel will play the role of discrete and continuous 
times, respectively. 



The factorization problem Given an element g £ G, in the big cell, and a set of deformation parameters 
s — {sa)aeSi t — itja)aes,jetf4seai ^+1 = {Oj Ij 2, • ■ • } N_i = {1, 2, ■ • ■ }, we consider the factorization problem 

Wo-g-Wg^ = S{s,t)-^ ■ S{s,t), S eG- and S €G+, (6) 

We will confine our analysis to the zero charge sector 

\s\ ■='^Sa = 0, 

aecS 

and consider small enough values of the continuous times. Observe that normally, but not always, the iofc times 
are disregarded. The reason is the triviality of factorization associated with this deformations. In fact, if we have 
a solution of the factorization problem for <ofe = 0, with factors Sq and Sq, then the factors corresponding to the 
factorization with arbitrary tok are S = exp{Y,k=i tokEkk)So exp(- X^f^i hkEkk) and S = exp(^f^^ hkEkk)SQ. 
The reason to consider them here is due to the reductions we will study later. 

At this point we discuss some relevant subalgebras and subgroups which will play an important role hereafter. 
Firstly, we notice that an operator A — J2jez ^ji^)^-'' commutes with A if and only if the coefficients Aj do not 

depend on n. Thus, the centralizer of A is 3a G g : [A, A] = 0} = | J2jez ^j^'' ^ ^ MAr(C)|. Observe 

that 3a C is a Lie subalgebra as now A commutes with the matrix coefficients of the Laurent expansions. 
Another interpretation is that we have block bi-infinite Toeplitz or Laurent operators [5] . 

A particular Abelian subalgebra f) of 3a is given by the centralizer of (^{A, Ekk}k=iJ i-'^i f) •= G g : 
[A, A] = [A,Ekk] = 0,fc = l,...,Af} = (Ejez^J^^^J ^ diag(A^, C)| where diag(iV, C) is the subalgebra of 
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diagonal matrices of Mn{C). Thus, [) is the set of Laurent series in A with diagonal n-independent coefficients. 
There are two important subgroups: G_ n 3a = {f + ciA^^ + C2A^^ + • • • , c j G M]y{C)} and G+ Ci — 
{cq + CiA + C2A^ + • • • , Co G GL{N, C), cj G Mjv(C), j ^ f }. Finally, we have the corresponding Abelian Lie 
subgroups H := G nt) and H± := G± fl () and Wq, Wq takes values in H. 

We shall denote by n e g the multiplication operator by the sequence {niN}nez', i-e. 

7i{X(n)}„gz = {nX{n)}nez- (7) 

Observe that [A, n] = A and that for any X G 2 we have X = ^^gz -^^yf^'A-', AT^j G MAr(C). This expansion 

follows from the assumption that Xj{n) = X^j + A"ijn + • • • . The set of operators commuting with A, n and 
Ekk, fc = 1, . . . , iV is given by {A G g : [A, A] - [A, n] - [A, Sfe,-] - 0, A: - 1, . . . , TV} = diag(A^, C). 

2 Lax and Zakharov— Shabat equations 

2.1 Dressing procedure. Lax and C operators 

We now introduce important elements for the sequel of this paper 
Definition 1. We define the dressing operators W,W as follows 

W:^S-Wo, W:=S-Wo, (8) 

In terms of these dressing operators the factorization problem(|6l) in G reads 

W-g = W (9) 

Observe that the expansions of the factors S, S 

S' = Dat + (pi{n)A^^ + ip2{n)K^'^ H G 

S = (foin) + (fi{n)A + (p2{n)k^ H G G+. 

Sometimes we will use the notation 

/3 :=</'!, e'^:=(po. 

We have the following expressions 



(10) 



N 

W ={^N + <^i(n) A-1 + (^2 (n) A-2 + ...). ^ Euk^'^ exp ( ^ t^k^^ 

k=l j=0 
N oo 

W = ((^o(n) + (^i(n)A + <f2{n)A^ + •••)•( XI ^kkA-"'^ exp ( ^ tf^A 



(11) 



k=l 3 = 1 



Other important objects are 
Definition 2. The Lax operators L, L, Cki, Cki,Cki,Cki S are given by 

L:=W-A-W-\ L:=W-A-W-\ (12) 

Cki -.^ W ■ Eki ■ W-\ Cm ■■= W ■ Eki ■ (13) 

Cm --^ S ■ Em ■ Cm '-^ S ■ Em ■ S~\ (14) 

Notice that in the above definitions of L,L,Ckk and Ckk — as Wq, Wq G H — we may replace the dressing 
operators W and by S* and S, respectively. A straightforward calculations yields 

Proposition 1. 1. The following relations holds 

oc 

Cm = L'"-'' exp(X(t,fe - t,i)L=))CM, 

OO 

Cm = L-^-^+''eMj2(t^k ' hT)^~'))^M. 
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2. The Lax operators have the following expansions 

i = A + + M2(n)A"^ H , Z"^ = uo(n)A"^ + ui(n) + U2(n)A H , 

Cki = Eki + Cfeia(")A"^ + Cfci,2(n)A-2 + . . . , Qm ^ CkiA^) + Cki,iin)A + Cki^n)^^ + • • • ■ 

3. These operators fulfill 

N N 

DAr = ^Cfcfc, iM = Y.^kk, (16) 

fc=l fc=l 

CkiCk'i' = Sik'Cki' , CkiL — LCki, ^^^^ 
CkiCk'i' = Sik'Cki', CkiL — LCki, 

2.2 Lax and Zakharov— Shabat equations 

In this section we will use the factorization problem ([9]) to derive two sets of equations: Lax equations and 
Zakharov-Shabat equations, and we will show they all are equivalent. Let us first introduce some convenient 
notation 

Definition 3. 1. 

d 



din. := 



' dt - ' 

2. The zero-charge shifts Tk for K = (a, h) are defined as follows 
and all the others discrete variables remain unchanged. 

3. 

QK TkWo ■ Wo'\ QK ■■= TkWo ■ K = (a, b). 

I 

n.a-.^wd.aW-^, nja:^weajW-\ Uk-.^WqkW-^ Uk-.^WqkW-K 



5. 



Notice that if 



7r„ := 



iOK ■■= {Uk ■ Uj,')- ■ Uk = {Uk ■ Uj,^)+ -Uk^G, {UkU],^)± e G±. 



Ekk, a^keS, 7f •= J^' 

0, a G S, ° 1 Ekk, a e S and a = fc for some /c e S. 



(18) 



(19) 



we can write 



gx = Djv + 7ra(A-DAr)+7rfc(A-i-Djv), = Djv + tt JA^^ - Dat) + 7f6(A - Dw), K = {a,b). 

Observe that all the shift operators preserve the zero charge sector and form a commutative group 

TkTk' = Tk,Tk, (20) 

T(a.b)T{b,a) = id, (21) 
satisfying the following cohomological relations 

T(a,b)T(b,c)T(c,a) = id. (22) 
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Proposition 2. The relations ([20| -([22ll are equivalent to 

T(a.b)T[b,c) = T(^b,c)T(a.b) = ^(a,c) (23) 

Proof. See Appendix B. □ 

Also notice that Bjk = {CkkL^)+, B^f. = {CkkL^^)- and that (HH]) and IH]) gives 

cjif = TTaA + OK + a^A^^ (24) 

for some matrix sequences axin) and aK{n). 

The factorization problem ([5]) implies that the partial differential equations 

djaW ■ = d,aS -S-^+S- d,a ' S'^ = O.aS • ^"^ + 5 • • = ■ (25) 

and partial difference equations 

TkW ■ W-^ = TkS ■ qK ■ S-^ = TkS ■ qK ■ S-^ = TkW ■ W"^ (26) 

hold. 

From the previous proposition we derive the following linear systems for the dressing operators and Lax 
equations for the Lax operators, and its compatibility conditions 

Theorem 1. 1. The dressing operators are subject to 

d,aW = ■ W, d,aW = B^a ■ W, (27) 

TkW = ujk-W, TkW ^lok-W. (28) 

2. The Lax equations 

djaL = [Bja,L], djaL — [Bja,L], djaCkk — [Bja,Ckk], djaCkk = [Bja,Ckk], (29) 

TkL — ujk ■ L ■ ujj/, TkL — ujk ■ L ■ ujk^ TKCkk = ■ Ckk ■ ^k, TKCkk = ■ Ckk ■ ^k^ (30) 
are satisfied. 

3. The following Zakharov-Shabat equations hold 

djaBib — dibBja + [Bib, Bja] ~ 0, (31) 
TKBja — djaUJK ' l^x^ + ' ^ja " ^ K^ : (32) 

Tkujk' ■ i^K = Tk'UJk ■ i^K'- (33) 

Proof. 1. First, observe that (^5]) implies djaS ■ S^^ + IZja = djaS ■ S^^ + Ti-ja, and therefore djaS ■ S^^ = 
— {TZja — TZja)- e 0- and djaS ■ = {Tlja — Tija)+ e 0+ so that using again (^5)) we get 

djaW ■ W~' = -{Uja - Tija)- + 7^,a = = (7^,, - 7l, J+ + Uja = djaW ■ W-' (34) 

Equation implies 

TkW ■ W-^ = TkS ■ qK ■ = TkS ■ 5"^ • Uk = TkS ■ 5"^ • Uk = TkS ■ qK ■ 5"^ = TkW ■ W''^ (35) 

so that_(Tif5 • S-'^)-^ ■ (TkS ■ S-^) = Uk ■ Uk^ and we conclude TkS ■ S'^ = [Uk ■ W^^)- ^ and 
TkS ■ S^^ ~ [Uk ■ 1^k^)+ £ G+ which introduced back in gives 

TkW ■ W-' = TkS ■ qK ■ S-^ = [Uk -Uk')- ■1^k=ujk 

= [Uk ■ Uk^)+ ■ Uk = TkS ■ qK ■ S'^ = TkW ■ W-\ (36) 
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2. From the definition (|12p we get 



TkL = {TkW ■ W^) ■ L ■ {TkW ■ W-'^) 



djaL=[d,aW ■W-\L], 

TkL = {TkW ■ W'^) ■ L ■ {TkW ■ W-'^) 



TKCkk - {TkW ■ W"^) ■ Ckk ■ {TkW ■ W-^r\ TKCkk = {TkW ■ W-^) ■ Ckk ■ {TkW ■ W-^)-\ 

and using ^ and ^ we find ^ and respectively. 
3. Tfie compatibility of (EH) and ^ imply (PT|) - (P5)) 



Observe also that the trivial flows iofcj A; = 1, . . . , are immediately integrated and if Lq, Lq, Cki and C'ki 
are the Lax and C operators corresponding to tok — for arbitrary tko we only need to conjugate these operators 
with exp(X;^^i Ekktko)- 

The compatibility conditions pip - (|33p for operators Bja and luk formally imply the local existence of a 
matrix potential <^ such that Bja = dja£. ■ and ujk — Tk£, ■ C^^; here, the potential ^ is a map to the Lie 
group G depending on the variables {tja,Sa} i. Moreover, any operator generates a gauge transformation so 
that Bja — > djaS, ■ + C ■ Bja ■ and ujk TkS, ■ ojk ■ providing new solutions of (l3T |) - (f33|) . 

Proposition 3. The relations 



We have seen that the Lax equations and Zakharov-Shabat equations PT|l - ([55)l appear as con- 

sequence of the factorization problem. The compatibility conditions for the Lax equations are satisfied if the 
Zakharov-Shabat equations hold. It is a standard fact in the theory of Integrable Systems that by construction 
the Lax equations imply the Zakharov-Shabat equations and therefore the system is compatible. In [30] is 
proven this fact for the differential equations (not the difference nor difference-differential equations) involved 
in the multicomponent 2D Toda hierarchy, that is that ([29|l =4> (|3T|) . Here we give an extended proof in order 
to include the continuous-discrete and discrete-discrete cases. 

Proposition 4. Let {L^Ckk}k=i '^"■'^ {LTCkk}k=i C g &e two sets, composed each of them of commuting 
operators, consider functions TZja G 0, Uk G G of L, Cu, . . . , Ckk and TZja £ S, Uk G G of L, Ch, . . . , Ckk, and 
define Bja andujK according to (fT9|l . then the Lax equations (pQ]) and (|30p imply the Zakharov-Shabat equations 

dSIll-dSSl). 

Proof. See Appendix B. □ 
2.3 The multicomponent Toda equations 

Here we write down some of the nonlinear partial differential-difference equations appearing as a consequence 
of the factorization problem Q . From ([M)) and ([55)1 , taking into account that S £ G_ and S € G+ , we deduce 
the following 

Corollary 1. We have the expressions 



□ 




(37) 



and the compatibility conditions (j33p are equivalent. 



Proof. See Appendix B. 



□ 



Bi, 



TTaA+Ua + UaA \ 



(38) 



COK ■— TTaA + OK + OK A 



K^{a,b), 
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where the coefficients have the alternative expressions 

Ua f3{n)TTa - TTaPin + 1) = 



0, a G S 



ax DtV - TTa - TTfc + TKl3{n)lTa - TTaPin + 1) 



0, aeS, 
diaP{n), a e S, 

''e^'^'^(")-(DAr-^fc)-e-'^("), aeS, 
Dtv - TTb, a e S, 



(39) 



\TK/3(n)(DAr - TTb) - (Djv - 7rfc)/3(n) + Tr^, a G S. 

From (|39p we deduce the following set of nonlinear partial differential-difFerence equations 

f3{n)Ekk - EkkP{n + 1) = 9u.(e^(")) • e'^^"), 
9u./3(n)=c^(")i?fcfce-^("-i), 
T(u,b)P{n)Ekk - EkkPin + 1) + - E^k - tt, = c^e^.'')^^") -(D^ - ir,) ■ e-^f"), 



(40) 



These equations constitute what we call the multicomponent Toda equations. Observe that if we cross the two 
first equations we get 

dik' {dikie*^"^) ■ e-^(") ) = e*(") E,,,,,, e-*("-i) E^k - E^k 6*^"+^) Ek'k' e'-^^") 
which is the matrix extension of the 2D Toda equation, which appears for iV 1: 

didl{(t)in)) = e'l'in)-4'{n^l) _ ^Hn+l)~4'{n) 

If in the last equation we set 6 = Z G 5 we have 

A(jj)/3(n) = e^(S.f)*(") -E,, ■ e-^^^-D . 

which when considered simultaneously with the first gives 

A(Pj)(aifc(e'^(")) • e-^(") ) = e^(^'.o-*(") -Ek'k' ■ e'^^""!) Ekk - E^u e^('^'.o^("+i) -E^k' ■ e'^^") 

which is a Toda type equation. A completely discrete equation appears, for example, when crossing the two last 
equations, i.e. 

A(fc^^)(e^<'=•')*(") •(Dw-^fc)-e-^(") ) = T(fc,b) ( e^(f -Ek'k'-e-^^"-'^ )Ekk~Ekke^cy .g-^W . 

So forth and so on we may get a set of continuous-discrete set of Toda type equations. Finally, observe that 
when iV = 1 we only have the shift Tj-^^ which corresponds to a shift n — > n + 1. 

2.4 Block Toeplitz/Hankel reductions 

We now consider some reductions of the multicomponent 2D Toda hierarchy. In the first place we discuss an 
extension of the periodic reduction [30 and the bigraded reduction 7 to the multicomponent case, which we 
call Toeplitz/Hankel reduction. Finally we discuss an extension of the 1 dimensional reduction discussed in 
[30] . These reductions are relevant when we work with semi-infinite cases, as in the construction of families of 
bi-orthogonal and orthogonal matrix polynomials, to be published elsewhere. 
Given a set {ta]a£S C Z we seek for initial conditions g satisfying 

N N 

g-[Y. E,,A~'-^) = ( J2 EkkA'-) ■ 9- (41) 

k=l k=l 
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The relation (|4T|) gives the foUowing constraints over the Lax operators 

N N 
k=l k=l 

for any j € I.. To proceed further in the analysis of these reductions we define the sets 

S± := {a e S : ±4 > 0}, So := {a £ S : 4 = 0}, S± := {a G S : ±4 > 0}, Sq := {a G S : 

so that S = S+ U So U S-, and S = S+ U So U S_. 
Proposition 5. // (j42|l holds then we have 

1. The dressing operators are subject to 

N N 

"aeS+USoUS+ k=l aeS+USoUS+ k=l 

N N 

aeS_USoUS_ fc=l aeS_USoUS_ k=l 

for j > 0. 

2. The Lax operators are invariant: 



( E 5A,a)(i) = ( E a,c,„)(Z) = o, 

aGS+USoUS+ aeE+USoUS+ 

( E 9,\,^Ua)(L)={ E a,|,„|.„)(Z)=0, 



aeS_USoUS_ aeS-USoUS- 

where j > 0- 
Moreover, if 

E = 0; 

aG5 

then, 

1. The dressing operators fulfill 

N 

k=l 

N 

W{S1 + ii, . . . , SN + ^N, Sl + if, . . . , Sj^ + i^) = Wisi, . . . , SiV, Si, . . . , s^) E EkkA~^ 

k=l 

2. The Lax operators are periodic 

L{si + ^1, . . . , Sat + ^AT, Si + ^i, . . . , Sjv + = i(si, . . . , sat, sj, . . . , Sjy), 

Z(si + ^1, . . . , Sat + ^AT, Si + €l, . . . , S^ + Ifj) = Z(si, . . . , Sat, Si, . . . , Sjv)- 

To prove this we need the 
Lemma 1. // (|42|) holds then for j > we have 

N N 

E = E '^kkL^"' = E CkkL-'", 

aeS+UEoUE+ k=l k=l 

N N 

E ^^K^u = E Cfefe^"'''' = E CkkL'-^'. 

aeS-USoUS- k=l k=l 
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Proof. The projection on g+ of A — J2k=i ^kkL^''-' , j > 0, is X^aes+uSo ^J^a,a while the projection on 0_ of 

A — X^aLi ^kkL~^>'^ , i > 0, is J2aes+ ^j^o.,a- The first formula is just A = A^ + A-. The second formula 
follows in a similar way when j < 0. □ 

Now we proceed with 

Proof of Proposition O Equations ()43|) and (I44p follow from the previous lemma and Theorem [1] To deduce 
(gSI and pS)) we argue as follows. If 

aes 

the periodicity follows from the factorization problem 

N N N 

k=l k=l k=l 

by observing that 

AT 

Wo{si +ei,...,SN +iN) = Wo(si, . ..,Sn)J2 Ekk^''"' 

k=l 
N 

Woisi + ^1, . . . , + ^jv) = Woisi, ...,s^)Y, EkkA-^'', 

k=l 



and recalling the uniqueness property of the factorization problem we deduce the periodicity condition for the 
solutions 

S{si + £i, . . . ,sn + iN,Si + £i, . . . ,Sf^ + %) = S{si, . . . ,Sn,Si, . . . ,Sf^), 

S{si + £i, . . . ,Sn + iN,Si + ii, . . . ,Sff + iff) = S{si, . . . ,SN,Si, . . . ,Sfi), 

which imply and ([46]). □ 

Now we justify the name of this reduction. If we write g = X^jez ■9i("')A"'i and think of it as an element in 
Mw(Mz(C)), i.e. g = Y.Zm=i 9kik2Ekik2 and gk,k2 = Ejgz Sj.fcifea then (gT]) gives 

9jMk2in) = 9j+iu^+iif,^Mk2{'n- - 4i)- (48) 

If ^fci + ^A:2 ~ ^' then gj.kik2 is a |-periodic function in n. If this period is 1, we get that gkik2 is a bi-infinite 
Toeplitz or Laurent matrix . We will see that in the general case we are dealing with block Toeplitz 5 and 
block Hankel [31] bi-infinite matrices. 

Definition 4. Given a block matrix = made up with p x q-blocks fii.j we say that fl is a block 

Toeplitz matrix i/ili+ij+i = fli^ and a block Hankel matrix i/f2i+ij_i — ^lij- 

Proposition 6. The condition ()48|) implies for g^-^^^ that 

• For Ilkiik2 > is a \tk-^ \ x Kfcjl"^^''''^ bi-infinite Hankel matrix. 

• For ik-^ik2 < is a \£ki \ x Kfc2 1"^^'-"^^ bi-infinite Toeplitz matrix. 

• For Iki = with i^.^ ^ we have a diagonal band structure being its width, and for (,^2 ~ with 
(■ki ^ a \(-ki \ X Kfcil block bi-infinite matrix. 

Proof. See Appendix B □ 
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The Toeplitz/Hankel block structure appears not only in the structure of gkik2 but also in the structure of 
g itself, thought as an element in Mz{M]\[{C)), for example if one takes £k = —^k — 1, k — 1, . . . ,N we get a 
N X iV-block bi-infinite Toeplitz matrix, while for = £k = ^ = 1, . . . , N we get a x A^-block bi- infinite 
Hankel matrix. 

Notice that for the particular case ik = ^ = 1, . . . , N , we have that g is a block Hankel bi-infinite matrix 
and 

N N N N 

g EkkA^'" = EkkA'^g, g Y ^kkA'" = ^ EkkA'^'g- 

k=l k=l k=l k=l 

From these two equivalent conditions on g we conclude for g^ the following constraint 

N N 
k=l k=l 

i.e. g^ is a block Toeplitz bi-infinite matrix and the corresponding solution to g^ of periodic type with bared 
and non bared periods equal to each other for each component: £^ = —ik, k — 1, . . . ,N. 
In the one component case we get the condition 

i^i=Z~^i. (49) 

li £i + il = we may choose £i — £ E \s4 and £i — —£ and the constraint for g is gA^ = A^g which leads 
— L^, i.e., the £-th periodic reduction of the one component 2D Toda hierarchy [30]. When £i,£i > are 
two nonnegative integers this constraint (|49p gives the the reduction of the one component 2D Toda hierarchy 
suitable to be extended with additional flows as described in [7J, named there as bigraded. This is why we refer 
to this reduction when all £a are positive as multigraded reduction. Notice that this multigraded constraint over 
g is never of periodic type and S = S+ and S = S+. 

ID reduction and generalizations Given a set of nonnegative integers {£a}a£S we request (7 in © the 
following constraint 

N N 

g. (5]£;,,(A^* + A-^*)) = (^i?,,(A^'= + A-^^)) ■ g. 

k=l k=l 

Now, as + z^^ = (z + z^^y + ajj^2{z + z^^y^^ + ■ ■ ■ + aj,o, for some aj^i eZ we have 

N N 



fc=l k=l 
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and therefore 

N N 



YCkk{L^'''+L-='^)=Y^kk{L='-^+L-^'-^) 

k=l fe=l 



is fulfilled for any j ^ 0. 

From here we conclude that 



N N N 



k=l k=l k=l 

and therefore we deduce the invariance 

N N 



J2^d,e,,k + = E(^A,fc + d^i,,k)L = 0. (50) 



k=l k=l 



In the one component case if we choose = ^ j = 1 we get the invariance under dji + dji, j > 0. This is the 1 
dimensional reduction as discussed for example in [30j . 
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It must be stressed here that being the same invariance conditions (|50p for this reduction and the previous 
multigraded reduction, the conditions are different for g and therefore for the class of solutions considered in the 
2D Toda hierarchy. In fact the Ueno-Takasaki ID-reduction has soliton solutions, which appear as a particular 
class of the general soliton solutions of his 2D Toda hierarchy. However this Ueno-Takasaki's family of soliton 
solutions of 2D Toda does not admit the bigraded type condition. On the other hand, the condition L = L^^ 
appears as a string equation in 2D Toda leading to solutions of the 1-matrix models, see for example [25]. 

3 Orlov— Schulman operators, undressing, and string equations 
3.1 Introducing the Orlov— Schulman operator 

Given solutions W, W of the factorization problem ([9]) we introduce the Orlov-Schulman operators [23] for the 
multicomponcnt 2D Toda hierarchy 

Definition 5. The Orlov-Schulman operators are defined as follows 

M:=WnW^^, M:=WnW-^, (51) 

Proposition 7. • The Orlov-Schulman operators satisfy the following commutation relations 

[L,M]^L, [M,Ckk]=0, [L,M]=L, [M,Ckk]=0, (52) 
• The following relations hold 



N oo 

k=i j=i 

N oo 

k=i j=i 



(53) 



Proof. See Appendix B. □ 
Given the initial condition g £ G m the factorization problem ([9]) we write 

N N 

gAEkk9~^ ^ ^ Pk,U'{'"-^^)Eii', gnEkk9~^ ^ ^ qk,w(n,k)Eii> 

IM = 1 l,l'=l 

and define 

N N 

Pk-.^ Pk,ii'{M,L)Cw, Qk^ J2 qk,u'{M,L)Cw, (54) 

= 1 l,l' = l 

SO that 

[Pk,Qk']--Skk'Pk- (55) 
Then, as Wg = W, we get, in the language of [26 , the string equations 

N N 

Pk,u'{M,L)Cw ■=LCkk, ^ -^^(^fefe- (56) 

i,i'=i i,i'=i 

3.2 Undressing Lax equations for the Lax and Orlov— Schulman operators 

The Orlov-Schulman operators M — WnW'^^, M = WnW^^ satisfy 

TkM = {TkW ■ w-'^)m{TkW ■ w-'^)-\ TkM = [TkW ■ w-^)m{TkW ■ w-^y^ 
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and the factorization problem ([9]) holds then the results of Theorem [T] imply the following Lax equations for 
Orlov-Schulman operators 

djaM = [B,a , M] , djaM = [B,a , M] , 

TkM = ujkM uj^^ , TkM ^ ijJkMuj^^ . 

We now prove the local equivalence between the factorization problem and the Lax equations. 
Theorem 2. Let us suppose that: 

1. The operators L,L,Ckk,Ckk,M and M satisfy the conditions 

L = A + ui(n) + M2(n)A"^ H , = ?2o(n)A"^ + ui{n) + il2in)A H , 

Ckk = Ekk + Ckk,i(n)A-^ + Cfefe,2(n)A-2 + ... ^ Ckk = Ckk,o{n) + Ckk,i{n)A + Ckk,2{n)A^ + ■■■ 

N oc N oo 

M = ■■■ + M_iA-i +n'Y^Ckk{sk + ^jtjkL^), M = • • • + Mi A + n + ^ Ckk{sk + Y.^*fkL~- 

k=l j=l k=l j=l 

with k — 1, . . . , TV, uo(n) G GL(7V, C), and fulfill the equations 

N 

Dw = ^ Ckk, CkkCii = SkiCkk, CkkL — LCkk, CkkM — MCkk, LM — ML, 

k=l 

N 

— ^ Ckk, CkkCii = SkiCkk, CkkL — LCkk, CkkM — MCkk, LM = ML. 

k=l 

2. Given operators B and lo as in (I19p . the Lax equations (|29p . (|30p and (I57p hold. 

Then, there exists operators S G G_ and S G G+ such that for W — SWq and W — SWq we may write 

L = WAW^^, M = WnW^, Ckk = WEkkW-^, 

L = WAW~^, M = WnW-\ Ckk = WEkkW-\ 

so that W and W solve the factorization problem ^ for some constant operator g. 

Notice that the set of constrains ([55]) and (|59p are preserved by the Lax equations. 
Proof. Observe thatwe need to find is the representation 

L^SAS-\ M^SnS~\ Ckk^SEkkS-\ 

L^SAS~\ M = SfiS-\ Ckk^SEkkS-\ 

with and p, as in ([55)1 . We first undress the Lax operators L, L^^. We look for 

5' = Dat + (p[{n)A-^ + (p2{n)A-^ + ■■■ , f']- ^ Mn{C), 

^' = (^;,(n)(Djv + <^'iHA + ^2 WA' + •■•), (^J, : Z-^GL(7V,C), : Z ^ Mn{C), J>0 

such that 

{A + UI+ M2A-1 H )(DAr + tp[A-^ + tp^A-'^ H ) = (D^ + tp[A"^ + tp'^A^'^ H )A, 

(■OoA^^ +Mi +U2A2 H )<^o(lljv + (p'lA + i^JjA^ H ) = Lp'oi^N + ip'^A + if'^A'^ H )A"^ 

Therefore, if we define T± [/] := X^jlo /(" ^ j)' have 

ip[ = ci + T+[ui]. 

Once is obtained we fix our attention in the second equation to get 

(f'^in) = C2 + I+[U2(") + Uitp[]. 
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So forth and so on we get all the coefficients Lp'^ up to integration constants Cj. 
Now we analyze (pT|) . which we write as follows 

(-SoA^^ + t;i + V2t^ ^ )(DAr + <^iA + i^'^t^ ^ ) + .^'i + (p'^lS. ^ — 

with Vjin) := i^o(n)~-'^Mj(n)i^o(^* + J ^ !)• From this we deduce that ?;o = or 'ito(?i) = (^Q(n)i^Q(n — l)^"'^- 
Denoting := logi^g S^* logwo = (1 — A^^)(0') and therefore ^'(n) = exp (cq + l^-[loguo]) where cq is a 
constant matrix. As uo = we have 

(A"i +vi+ usA^ + . . . )(D^ + ip[A + (p'^t^ + •••)= A"^ + (^'i + (^^A + • • • 

which transmutes into (|60p once we replace Uj by 5j, i^a^- by i^^-, j — 1,2, .. . and A by A^^. Thus, all the 
coefficients (p'^ are expressed in terms of Vj . 
Now, we proceed to undress C^fc and C-k,^^ 

These operators commute with A, satisfy 



C'kk - Ekk + C'kk,!^ ^ + Ckk,2^ ^ H 

C'fcfc = C'fcfc.o + C'fefc^iA + Ckk^2^'^ + • • • 1 



and provide us with two different resolutions of the identity, 

N N 



Hjv — ^ Ckk, CkkC'ii — SkiC'kk, Dat — ^ C'kk, CkkC'ii — SkiC'k- 

fc=l k=l 

In fact, it is easy to show that there exists operators Q — ^n + Qi + • • • G G- Dja and Q = Qo + Qi A + ■ • • G 
G+ njA such that C^^, = QEkkQ~^ and (7^,;, = QEkkQ^^ ■ Thus, to undress L,L,Ckk,Ckk, fc = 1, . . . , iV, we 
just take S* = 5" • Q, = 5' • Q. 

With these operators at hand we proceed to undress M and M 

S-^MS ^ a + jjL, a:^ S-^MS -n, S^^MS ^ a + jl, d:^ S'^MS -n, 

but 

[A, S-^'IS] = A, [Ekk, S-^MS] = 0, [A, S-^MS] = A, [i;fcfc, ^-^A/S] = 

and [A, /i] = A and [A, p] — A, so that 

[A, a] = 0, [Ekk, a] = 0, [A, d] = 0, [-Bfcfc, a] =0 ^ a, a G f). 

Now, recalling that A4 = n + and n + Ag+ we write a — S^^nS — n + g_ and d — S^^nS — n + g+A 

so that a = aiA^^ + a2A^^ + ■ ■ ■ and d — diA + d^lS? + • • • with a;, di G diag(A^, C) for all i G N. We define 
7 - Ej>i T^"'' 7 00 + Ej>i ^A^ where 0o G diag(iV, C), and find that 

e''' n e'^ ~ n + [7, n] = n + a = S^^MS, e^ n e""*" = n + [j,n] — n + d — S'^MS 

which allows us to write 

e^ WquWo^ e"^ = e'^ ne"'^ +iy = S'^MS + e^ W^nW^^ e"^ = e^ ne"^ = MS + P. 

Therefore, if we replace 5 — > S'e''' and 5* ^ S'e''', we get the desired result. 
From the evolution equation we get that 

{Sy^ ■ {B,a - djaS • {S)-')S, A,a := {S)-' ■ {B,a - d,aS ■ {S)-')S, 

_ - - ybl) 

Pk-^Tk{S) ujkS, pk-=Tk{S) lurS, 
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commute with A and all the Ekk, k = 1, . . . , iV; i.e., they are n-independent and diagonal. We may deduce that 

[Ajn,^j] = jOja ^ [■Aja,n\ = j6ja [A-ja ^ dja,n\ = 0, 



(63) 



and 



which imply 



Thus, 



[Ajn,lj] = -jOja [Aja,n] = -jOja ^ [Aja " Oja,n] = 0, 

ujkMw],' = {TkS ■ S-^)M{TkS ■ + {TKS){^a - 7Tb){TKS)-\ 

QkMuj^^ = {TkS ■ S-^)M{TkS ■ S'^)-^ ~ {TKS){^a - ^b){TKS)-\ 

[pk, = (T^a - T^b)Pk ^ [Pk,n] = (vTa - Trb)pk ^ n] = 0, 



(64) 



Aja - Oja, Aja - Oja, Pkix , Pkix ^ diag(iV, C). (65) 



As the Lax equations (|29|) and (|30|) are satisfied by Proposition |4] we know that Bja and satisfy the 
compatibiUty conditions ((3T|l - ((37|) . However, we see from (|62|) that {Aja,PK} and {A^a, } are gauge transforms 
of Bja, ujk and thereby do have zero curvature. Therefore, we conclude the local existence of potentials f and ^ 
such that 

A,a=djai-C\ Aja=djabt\ PK=TK^-r\ PK = Tk^ ' ■ (66) 

These potentials are determined up to right multiplication C — >■ ^ • /i, f — > ^ • /i, where h,h £ G are constant 
operators independent of tja, Sa- Up to this freedom we may take the potentials ^, ^ G H.Now, recalling (fTO|) we 
get 

which together with ((^ imply — ^'ja G 0-, ^ja — ^ja G 0+, PkQ^^ G G_, pkQk^ G 6*+; but these operators 
belong to diag(A^, C). Thus, from the two first relations we conclude Aja — Oja, pK — qii and ^ = Wq while the 
two second imply that if ^ = e*^" -Wo then 0o G diag(A^, C). 
Therefore, we may write 

A,a = d,aWQ-W^\ A,a = d.amo) ■ {e't"> Wo)-\ pk=TkWo-W^\ =T;f(e^°W^o)•(e^°VKo)"^ 
We make the replacement S ^ S e"^° to get 

Bja = djaW ■ = djaS ■ S'^ + S0jaS-^ = djaS ■ + SSjaS-^ = djaW ■ 

UJK = {TkW) ■ W-^ = {TKS)quS-^ = {TKS)qKS-^ = {TkW) ■ W-\ 

In terms of 5' = ■ W the previous equations can be written as djag — and Txg = g- Thus, we finally find 
Wg = W where 17 is a constant operator in G. □ 

A further result regarding the operators Gki introduced in Definition [5] and characterized in Proposition [T] 
that will be needed later is given now. 

Proposition 8. Given operators L, L, M, M ,Ckk and Gkk as in Theorem\^ then: if we find operators Cki of 
the form 

Cki=L"'~'' e('^'=-'^')^'(Sfe, +0_), 

such that 

[Cki,L] = [Cki,M] = 0, Ck'k'Cki — Skk'Ck'i, CkiGk'k' = ^ki'Ckk', 

then the undressing operator W of Theorem\^ satisfies Cki — W EkiW~^ . 

Proof. See Appendix B. □ 
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3.3 String equations, factorization problem and Lax equations 

We will show here that the string equations (|56p for the Lax and Orlov-Schulman operators do indeed imply 
the factorization ^ and also the Lax equations In fact, only one of these implications is needed as 

the other one will follow from the results described previously. However, we show that these two facts can be 
derived directly from the string equations, showing the importance of this formulation of integrable systems. 

Theorem 3. Let L, M,Ckki L, AI, Ckk, ^ — l,...,iV, be operators as in Theorem \^ and operators Cki, k,l = 
1, . . . ,N, be as in Proposition\^ Let us suppose that we have operators Pk,Qk os in (I54p and that the string 
equations (|56p hold. Then, 

1. We can choose the operators W and W of Theorem\^ such that the factorization ^ holds for some constant 
operator g d G. 

2. The Lax equations (|29 p - (|30p are fulfilled. 



Proof. From the first part of the proof of Theorem [2] (not considering the Lax equations) and Proposition [8] we 
know that there are undressing operators W — SWq and W = SWq, S £ G_ and S £ G+. Let us introduce 
some convenient notation 



OK {TrW ■ W-^)-^TkW ■ W-^ a% := W'^grW, 



(67) 



and observe that if we define 



( W-^ ■ W (68) 

we have 

D%=d,a(:-C\ a%^TKC-C-'- (69) 

The string equations ((56)) read 



N 

Pk = WPj^W"^ = WEkkAW~\ P^ = PkM'{n,K)Eiv 

i,i'=i 

N 

Qk = WQlW-^ = WEkknW~\ Ql = ^ qkw {n, A)Eu' 

Ll' = l 



(70) 



which in turn imply 

djaPk = [djaW ■ W~\Pk] = [djaW ■ W-\ Pk], 
djaQk = [djaW ■ W-\Qk] = [d,aW ■ W~\Qk], 

TxPk = {TkW ■ W-^)Pk{TKW ■ W-^r^ = {TkW ■ W-^)Pk{TKW ■ W~^)-\ 
TKQk = {TkW ■ W-^)Qk{TKW ■ W'^Y^ = {TkW ■ W-')Qk{TKW • W^Y^- 

Hence, recalling Pk — LCkk and Qk ~ MCkk we conclude that 

[Dja,L] = [Dja,M] = [D,a, Ckk] = 0, [^K, L] ^ [aK,M] = [aK, Ckk] = (71) 

Thus, we deduce that D'j^,(j'^ e diag(A^, C), and therefore Dja G 0+ and aK G G+. With these preliminaries 
let us start proving the two statements in the theorem 

1. Given the representation in terms of ( as in ([55)1 we deduce that we can write C = £, ■ for some 
^ e diag(Af, C) and some (s, t)- independent operator g e G. Thus, = djaS, ■ and cr^ = TkS, ■ ■ 
But, recalling the definition (|S5)) we get WC, = W ^ SW^S, — Wg. Observe that [Wo,C] = and replace 
S ^ S ■ £, to get the factorization problem. 
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2. From the definition (|67p we get 

TkW ■ = TkS ■ ■ Uk = TkS ■ ■ Uk ■ ctr = TkW ■ W'^ ■ ax- 
Reasoning as in the proof of Theorem [5] and recalhng that Dja G 0+, cfk G G+ and [W/f , (Jk\ = we have 



SO that, according to (fTO|l . 

a^aM^ • W^"' = - a^al^ • W^'^ + D,a, TkW ■ W'^ =ujk = TkW ■ W-^ ■ gk- 

Therefore, we immediately get the following Lax equations 

djaL = [d,aW ■ W-\L] = [Bja,Ll Tk L = [TkW ■ W-')L{TkW ■ W-^)-' = ukLlu],' 

djaM = [djaW ■ W~\M] = [Bja,M], TkM = {TkW ■ W-^)M{TkW ■ W'^y^ = ujkMuj^^ 
djaCkk = [djaW ■ W~\Ckk] = [Bja,C'kk], TkCuu = {TkW ■ W-')Ckk{TKW ■ W-^' - ujKCkkC^K'- 

Now, as djaW ■ W~^ ~ Bja — Dja and TkW ■ W^^ = ujk ■ <^k with Dja and (Jk commuting with any 
function of L, M and Ckk we get the remaining Lax equations 



djaL = [djaW ■ W~\L] = [B,a,L\, Tk L = {TkW ■ W-')L{TkW • = ukLu]^' 

djaM = [djaW ■ W-\M] = [Bja, M] , TkM = {TkW ■ W-^)M{TkW ■ W'^^ ^ lukMlo]^^ 
djaCkk = [djaW ■ W-\Ckk] = [Bja,C'kk], TKCkk = {TkW ■ W-')Ckk{TKW ■ W'^ - coKCkkiOj,'. 

□ 

The above result might be slightly generalized by considering string equations of the form 

N N N N 

9fe,H'(Af, L)Cu>, 

i,i'=i i,i'=i i,i'=i i,i'=i 

where we assume that 

AT N 

i,i'=i i,i'=i 

belong to the adjoint orbit of Ekk-A, Ekkn, k — 1, . . . , TV; i.e., there exists g & G such that 

P°=g-EkkA-g-\ Qt=g-Ekkn-r'- 

For that aim the proof needs to be modified only in the definition of D'ja g^^ D^jag^^ and ct^l' ~* 9^^'^k9 
9 ~^ 9 ' 9- Observe that elements in can be constructed in terms of operators Ck,V and Q, such that: 

N 

Y Ck = Dw, CkCk' = Skk'Ck, [Ck,r] = [Cfe, Q] = 0, [V, Q] = V. 

k=\ 

What we do not know yet is if the orbit G is characterized precisely by this properties. However, if we request 
the following properties: Ck — Ekk G 0- , 7^ — A e g- A, Q — n G 0- , one could prove, following similar arguments 
as in the proof of Theorem [2l that these elements lay in . This implies an alternative formulation of string 
equations ((56|) 

Y^k,iv{L,M)Cu' = Ckk, ^7'ir(i,M)Gr = Z, Y Qu'{L,M)Cu' = M. 

1,1' 1,1' 1,1' 
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3.4 Additional symmetries and string equations 
3.4.1 Additional symmetries 

Suppose that the operator g in ([5]) depends on an additional, or external, parameter 6, which might belong to 
C or to Z. Now, we will describe the induced dependence on the elements defining the multicomponent Toda 
hierarchy. We shall denote by dt = djdb when £ S C is a continuous an by Tf, the corresponding shift £ ^ * + 1, 
when * € Z is an integer. In this case, we shall replace ^ by the equivalent factorization problem W ■ h = W -h, 
with 

g = h-h-\ (72) 

Observe that for 6 e C we may write, 

di,W ■ W'^ + W{dth ■ h~^)W^^ = d,S ■ + W{d,h ■ h-^)W-^ 

(73) 

while for b ^'Lwe have 
T,W ■ W-^ ■ W ■ {T,h ■ h-^) ■ W-^ = T,S -S-^ -W ■ {T,h ■ h'^) ■ W'^ 

= TeS -S-^ -W ■ {Tth ■ h'^) ■ W-^ = TeW ■ W^^ ■ W ■ {Teh ■ h"^) ■ W^^ 

(74) 

Now we suppose that the dependence on 6 is given by the following equations 

N N 



d,h -h-^ =Fo= ^ Fw{n,A)Ew, d,h ■ h'^ = Fq ^ ^ Fw{n,A)Ew when 6 e C, 

i,i'=i i,i'=i 

N N 

T,h ■h-^=To= ^ii'{n^^)Eiv, T,h -h-^^To^ ^u'in,^)Ew,, when B e Z, 



(75) 



and define 



Ll' = l Ll' = l 



N N 



i,i'=i 

N N 



(76) 



^ .F/i.(M,^)Q;', ^ J^iF(M,Z)Q;-, when £ £ C. 

i,i'=i i,i'=i 



From we get 

d,W-W-^ ^dr,S-S-^ = -{F-F)-, dsW-W-^^dsS-S-^^iF~F)+, {F - F)± e e±, 

T,W -W-^ =T,S ■ S-^ = ■ f-^)-, TiW ■W^'^ =T,S ■ S-^ = {T {T ■ T-^)± G±. 

So that 

Proposition 9. Given a dependence on an additional parameter b according to (j72p and (j75p . introduce H :— 
F ~ F and TL :— J- ■ where F and F arc defined (j76p . then 

1. The dressing operators W and W satisfy the following equations 

dsW = -H^-W, dtW = H+-W, or TiW^H-'W, TM ^n+-W. 

2. The Lax and Orlov-Schulman operators are subject to 

d,L = [-H^,L], d,M = [-H-,M], d.Ckk = ~\H^,Ckkl 

d,L = [H+,L], d,M = [i?+,M], d.Ckk = [H+,Cuk], 

or (77) 

T,L = H--L-HZ\ T,M = n--M-HZ\ T.C'kk = H- ■ Ckk ■ 'HZ\ 

T[,L = H+ ■ L ■ TLj^^ , TfM = H+ ■ M ■ 'H_^^, T^Ckk = 'H+ ■ Ckk ■ 'H_^^. 
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3.4.2 String equations 



The factorization problem ^ depends decisively on the 'initial data' g. Now, we are going to see some conse- 
quences of the form of g and derive the so called string equations. Let us suppose, that given two operators 

N N 

i,i'=i i,i'=i 
we have the following constraint satisfied by g 

9F0 = Fog. (78) 

Then, if 

N N 

F{M,L) := Fu>{M,L)Civ. F{M,L) = ^ Fu'{M,L)Cu', 



Ll' = l LI' 



we have 



F{M,L) = F{M,L). (79) 



We refer to these type of equations as string equations, see for example ^25], and we have seen that they 
reflect properties like ([75)) of the initial condition g in © . Notice that the reduction of (|¥T|l is a particular case 

of ([78| with Fq := X^feLi ^fefe^"^'' ^^'^ -fo '■= SfcLi -^fcfc^^'°: which the Orlov-Schulman operator does not 
appear. This suggests an important family of diagonal string equations with 



N N 



Fo ■.= Y,FkkFoAn,A), Fo ■.= J2EkkFoMn,A)- (80) 



fe=i fc=i 



The equations (|56p are also a set of string equations; moreover, the invariance conditions under the additional 
flow ((77l) implies that iJ = or 7i = id so that we are lead to the string type equations of the form (|79|) . namely 



F{M, L) = F{M, L) or T{M, L) = T{M, L). (81) 



This also follows from 



dtg = {dih ■ h ^)g - g{dth ■ h ^) ^ Fog - gFo, 
T,g = {T,h ■ h-') ■ g ■ {T~h ■ h-^' - -^o • .9 • ^o'. 

Observe if we consider arbitrary forms of Fo, -Fo or jFo, J-q it will be same to deal with the description given 
here or the one obtained just setting h — id. However the situation is different if we consider the function Fo, 
Fo or To, To of diagonal type ((80)) . In this case, to set /i = id will generically imply to abandon the diagonal 
family for Fq. 



Appendix A: Congruences 

We will show here how to derive from the multi-component Toda hierarchy equations involving only the fields 
at each site n G Z; i.e. not mixing fields at different values of n, the sequence variable. This is particulary useful 
to show the role of the discrete multicomponcnt KP hierarchy in the multicomponcnt Toda hierarchy, which 
appears when we froze the bared continuous and discrete times. 
Firstly, we present a queue observation 

Proposition 10. Let us suppose that we have operators R, R € q such that 

RWo-'€S-, RWo-'e2+ (82) 

satisfying R ■ g — R. Then R — R ^ 0. 
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Proof. We have 



R^Rg = RW^^Wg = RW-^W => RW-^ = RW"^ 

therefore RW^^^S''^ = RW^^S-^ and recaUing and the fact that S e G- and S E G+ we conclude the 
statement. □ 

Next, and without proof (which consist in a systematic and sometimes elaborated application of the previous 
result) we show the appearance of some well known integrable hierarchies within the multicomponent Toda 
hierarchy. We firstly point out that continuous variables, and for each value of n, we have solutions of the 
TV-wave hierarchy and its modifications, moreover some discretizations of the modified iV-wave equations are 
proposed. These results are just a manifestation of the fact that if we froze the bared times we are just dealing 
with a discrete TV-component KP hierarchy in the spirit of [1]. Next, we recover within this context the the 
quadrilateral lattice equations. Finally, we present what we call the dispersive Whitham hierarchy in complete 
analogy to the one proposed in p7]-P5]. 

iV-resonant wave equations and its modification We introduce 

9:=5iiH hSijv, d:^dii-\ K^i^, 

in terms of which we have 

Theorem 4. The dressing operators satisfy the following equations 

djaW = Q,a{W), djaW^QjaiW), (83) 

where 

Q]k = Ujkjd^ + Ujk,j-id^~^ H h Uj7c,o, Qjk = '"jk.jd^ + Vjk,]-id^~^ H h v.jk,id 

with the coefficients Ujk,i,Vjk,i depending on dfipr,df(pr , respectively 



Ekk, i = 0, 

ViEkk - 21=0 '^]k,]-a<yj-a,i-a, « = 1, • ■ ■ , j, 



(paEkkO-.g, i = 0, 



" J2a=0'"Jk,j-a<Jj^a,i-a]EkkCrj]:iQ, i = 1, . . . , j - 1, 



and 



'J 



r=0 ^ ^ 

^J,^■■^i2(^)i9''^^-r). (85) 



r=0 ^ ' 

Observe that Uji — (3 and aj^ — e"^, and also that the first of the differential operators Qjk and Qjk are 
given by 

Qjk = Ekkd' + [/?, EkkW'^ + {[^2,Ekk] - jEkkdP - Ekkl3\)d'-'' + ■■■ + Ujk,o. 
Qfk = Ekk e'^ + {^A - Ekk e'^i^i + jde*)) B^-^ + ■■■ + v,k,id. 

Lemma 2. The only differential operators Q = Ujd^ + ■ ■ ■ + uq and Q — Vjd^ + ■ ■ ■ + vq such that 

Q{W) = 0, Q{W) = 

are 

Q = g - 0. 
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Proof. • Let us suppose that J2t=a Uid'W = but X^Lo Uid'-W = ( I]Lo("^+X]r=i Ui+r<yi+r,r)^' +Q-))Wq. 
Thus, Uj = Uj-i = • • • = Mo = 0. 

• Assume now that X]i=i vtd^W — and take into account that 

^^d'W = ( v,+ra,+,,r)A-' + Q+)Wo. 

1=1 1=1 r=0 

Thus, Vj = Vj-1 — ■ ■ ■ — vi = 0. 

□ 

From this lemma it foUows that 
Proposition 11. The Zakharov-Shabat conditions 

djkQii - diiQjk + [QihQjk] = 0, dfkQa - d^Qf^ + [Qa, Qf^] = 

hold. 

Proof. Just consider the compatibihty conditions of ((83|) together with Lemma [2] □ 

The site independent relations described in Theorem 3] constitute the A^-wave hierarchy, for the non bared 
flows, and its modification for the bared flows. These multicomponent equations contains may integrable systems 
[17], for = 1 we have the KP equation in nonbared variables and the modified KP equation in bared variables, 
for N — 2, the Davey-Stewartson equation in the t-variables and the Ishimori equation in the i-variables. For 
N — S we find the 3-resonant wave system (i- variables) and a modified version of it (t-variables). 

Proposition 12. The N-wave equations 

dik[f3, Eu] - du[f3, Ekk] + [[P, Eul [/3, Ekk]] = 0. 
and the modified N-wave equations 

dik{vi) - diiivk) + vid{vk) - Vkd{vi) = 0. 

with Vk '■— ipoEkk'fiQ^ are satisfied. 

Proof. The A^-wave equations appear as the compatibility of the Qik — Ekkdi + [/3, Ekk]- A "modified" A^-wave 
system [17 appears for when one considers the compatibility Qife = Vkd, k — 1, . . . , N . □ 

Discrete versions of the modified A^-wave equations For a fixed Z = 1, . . . , A^ let us introduce the 
following shift operator 



k=l,..^,N 
k^l 



and the operator 



XiA) ■.^J2^(kj){A)Ekk, Pi AeQ. 

k^l k^l 

Finally we also introduce, 

nk:^VEkkV-\ V ■.^En+Xie'^), 

and the difference operators 

\kj) n-kj) - 1, ^ ^= E ^f-{N-l). 

k^l 

Proposition 13. The dressing operators W and W satisfy 
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Conjugate nets and quadrilateral lattices We show now the role of conjugate nets and quadrilateral 
lattices as a part of the multicomponent Toda hierarchy. For this aim we first prove 

Proposition 14. If Ci G Mjv(C), i = 1,2 are such that eiEkk = £2Eii then 

ei{dik - l3Ekk)W = e^{TK - [Tk^Eu + Djv - Eu - 7r„))W, K = {I, a) 

eiidik - pEkk)W = e2{TK - {Tk^Eu + Djv - Eu - na)W, K = {I, a) 

ei e-"^ d.kW = 62 e-^-^A^ + ttJT^, K = {l, a) 

ei e-* d.-^W = €2 e-^-^Ax + na)W, K = ([, a). 
A particular consequence is 

^iiV2Ekk + 9ikVi - fiEkkfi) = £2{T(k,i)'P2Ekk + 7fei<^i(Djv - Eu - Ekk) 

+ Ell — {TkifiEkk + Djv — Ell - Ekk)'fi)- 
If we right multiply this relation by Em'rn' with m! ^ k and we take 

1. ei = Emm with m ^ k and €2 = 

2. ei = and €2 = Emm with m^k,l 

3. £1=0 and £2 = En 

4. ei = £2 = Ekk 
we find 

dlkPmm' - PmkPkrn' = 0, for 771, m k 

\k,i)Pmm.' - {T(kd)Pmk)Pkm' = 0, for m, m' ^ k, I, 

{T(k,i)Pmk)Pki +/3mi = 0, rriy^ k, I, 

T{k,l)(^lm' - {T(k,l)0lk)Pkm' =0, m' 7^ k, I, 

{Tik,i)Pik)Pki -1 = 0, 
dik log/3fem' - "^''^i^^^"' + \k,i)f3kk = 0, m' 7^ /, 

dik log/4; + ^{k,i)Pkk = 0. 

The dispersionfuU Toda— Whit ham hierarchy We fix Z e S and consider the shifts T(^a,i) for a e <S 
with a ^ I, and as we can not put a = I for a' G S, a' ^ I, 

Proposition 15. 1. For a', I, a' ^ I, there exists scalar operators 

^Jl = ^(la') + ^3l,0-lT(iy) H h Bjifl, 

aji = d'ii + aj/j_2a^r^ + • • • + aji^o 
where the coefficients Bji^i and aji^i are scalar polynomials in the T^i g^iy shifts or the du- derivatives of 
l3n, ip2,u V>j,n, respectively, for example -Bj/j-i = Pn - T^i^^^'^Pu such that 

djiiEuW) = ^jiiEnW) = ajiiEiiW), dji{EuW) = ^ji{EuW) = ajiiEnW), (86) 

2. For a ^ I there exists scalar operators 

where Bja,i are scalar polynomials in the T(i^a)-shifts ofPik, ^2,ik, ■ ■ ■ , 'Pj,ik when a = k and of(po,ik, ■ • ■ , <fj-i,ik 
for a = k, for example 



such that for a ^ I 



Ik 



Tlk,i)iPik) 

'fiO,lk 

^lk,i)^^o,iky 



a = fc e S, 
a = fc G S, 



dja {Ell W) = i^ja {Ell W) dja{EllW)=^^ja{EllW). 
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Appendix B: Proofs of Propositions 

• Proposition [2] Obviously is implied by (gOll-llSl)- It is also easy to conclude that and ([221) follow 
from (|23p . The non trivial part of the proposition is to prove that implies ([^0)1 : 

^(a,b)-'-(c.d) — J- (a.c)-^ (c,b)-l^ (cfi)-'- (b,d) — ^ (a .c) -L (c .b) -L (b ,d) (c .b) 

— J-(a,c)-'-(c,d)-L(c,b) — J-(c.d)-L(a,c)-L(c,b) — ^ (c .d) -L (a ,b) ■ 

• Proposition [3] We only need to show that (|37p implies ([55)1 as the reverse is evident. We proceed as in 
the proof of Proposition [2] 

{Tl^a,b)^{c,d))^{a,b) = {T(a,b){T{b,d)^(c,b))^{b,d))){'^( 

= {T(a.b){T(b,d)^(c,b))){T(a,c){T(c,b)^(b,d))^(^c,b))^{a,c) 
= {T(a.b){T(i,^4)UJ(cfi))){T(a,c)^(c,d))^(a,c) 
= {T{c,d)iT(a.c)^(c.b))){T(c,d)^(a,c))^(c.d) 
= {T{c.d){T(a.c)^(c.b))^(a.c))^(c,d) 
= {T(c.d)^(a,b))^(c,d)- 

• Proposition [4] We do not prove the differential case and refer the reader to [30]. Therefore we proceed 
to the remaining cases involving discrete times. 

1. We start by proving ([5^ . First, from the definition (fTO|) we deduce that 

d.aUJK ■ LOk^ = djaiUK ' U^^)- ■ {Uk ■ Uk'T- + [Uk ■ UK^)-djaKK ■ U^^ {Uk ■ Uk')Z^ 

= djaiUK ■ Uk^)- ■ {Uk ■ Uk')-' + i^i< ' ^x')-(-Bja - UkB^Mk^Wk ■ Uk')-' 
= d,a{UK ■ Uk^)- ■ {Uk ■ Uk')-' + {^k ■ Uj,^)-B,a{UK ■ Uk')-' - ^KB^aCOj,' 
and similarly 

= dja{UK ■ Uk')+ ■ {l^K ■ Uk') + ' + {^K ■ U]^^) + Bja{UK ' Uk') + ' " ^K-BjaW^' 

so that 

djaLOK ■ LOj^^ + UJKB.aUJK' = ^..(Wk • Uj^^)- ' {Uk ' W^')!' + {Uk ' Uj,')-B,a{UK ' Uk')-' 

= dja{UK ■ Uk')+ ■ ("a- • Uk') + ' + {^K ■ Uj,') + B,a{UK ' U]i%' ■ 

Now, using p3p and the commuting character of the Lax operators we get 

TKBja — {Uk ■ Uj(^)-Tlja{UK ■ Uj^^)_^ — TK{TZja — iZja)- 

= {Uk ■ Uj,')+nja{UK ■ Uk')+' + TK{nja - n^a)- 

and we deduce for / :— dja^K • + ^KBja^]^^ — TKBja the following expressions 

/ = dja{UK ■ Uk')- ■ {^K ■ Uk^)-^ - {Uk ■ UK^)-{Tlja - 1lja)-{UK ■ Uk')-' + TK{n,a - Hja)- 
= dja{UK ■ Uk^)+ ■ {Uk ■ UK')t + {Uk ■ Z^^')+(7eja - ^,a)+(WK • Uk') + ' - TK{nja ~ nja)+. 

which hold only if / = 0, as desired. 

2. Let us now prove From ^ and ^ we get 

Tkujk' = Tk{Uk'Uk')- ■ ^kUk'UJk' = Tk{Uk'U^})+ ■ ujkUk'I^k' 
or using (|19p again 

Tklok' ■ coK = Tk{Uk'Uk})- ■ {UkUj,^)-UkUk' = Tk{Uk'Uk')+ ' {UkUk')+UkUk' ■ 
Then, we deduce 

{Tkujk' ■ i^K)+ = {UkUk')+, {Tkujk' ■ (^k)- = {U^}U]^^)- 

Interchanging K ^ K' and recalling the commuting character of the Lax operators we get the desired 
result. 
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• Proposition [6] Let aij denote the elements of the bi-infinite matrix gkik2-, we now proceed to analyze the 
meaning of in different situations: 

— Block Hankel case Let us assume that iki^k2 ^ particular let us discuss the case where both 

integers are positive. If we start from the element aij the equation (|48p says that it is equal to some 
other element. To determine the this element in the matrix we observe that (j48p requires to move in 
the i-th row iki + ^k^ positions to the right and in the diagonal passing through that position go left 
^fej positions on this diagonal, i.e. go up tkx positions and to the left also positions. This gives us 
the block structure over off-diagonals as illustrated below. 




For negative integers , ^k^ ^ ^ have a similar discussion, replacing right motions in the row with 
left motions and up motions in the diagonal with down motions in the diagonal, and the same block 
Hankel structure appears. 

— Block Toeplitz case We now assume that tki^k^ < 0. Suppose that is positive. Then, when 
Ik^ > l^j.^! if we start from the element the equation (pS)) says that it is equal to some other 
element, say ai'ji. To determine the row i' and column j', we observe that (j48p tell us to advance in 
the i-th row tki — \^k2 1 positions to the right and in the diagonal passing through that position go left 
£ki positions on this diagonal, i.e. go up £k-i positions and to the left also ik-^ positions. So that we 
have 

fly = a,_^^^ = ai_£^^ 

For the case 4^ < |4J we use 5j,fcife2("-) = .^j+Kfc^ l-f^i ,fcife2 + 4i) so that we move |4J - 4i 
positions to the right and on the diagonal ik-^ positions down, which amounts to Ik-^ rows down and 
£ki columns right, i.e. 

and we get the same result, which immediately tell us about the block structure over diagonals as 
illustrated below. 




A similar discussion goes on for the case of negative iki and positive i^.^ . 
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— The case £ki — with £^.^ ^ gives gj,k^k-i(n) — gj+i-^.^.k^k^iin), which implies diagonal band structure, 
whether for i^^ = with 4i 7^ gives gj^hi^n) = 5^+^^^ ,fcifc2 + 4i), which describes a Ik^ x 4i 
block structure. Notice that these two cases can only exist for two or more components. 

• Proposition [7] Let us compute 

M = WnW-^ = SWquW^^S-^, (87) 
M = WnW~^ = SWonWo~^S-\ 

for this aim we must take into account that 



N 



fx := WonWo ^ = n + v, v = ^Ekk{sk +^jtjkA'), 

k=i j=i 



N 



(89) 



fi:=WonWQ ^=n + iy, v ^ -^Ekk{sk + ^jtfkJ^ 

k=i j=i 



Therefore, from (gT]) and ([551) we deduce that 

N oo 

k=l 3=1 
N oo 

M = SflS-' = SnS-' ~ J2 Ckk(sk +J2jtfkL-n- 

k=i j=i 

Finally, 

M := SnS-^ = (1 + /3(n)A-i + ip2{n)A-^ + ■■■ )n{l + (3{n)K-^ + ip2{n)K-'^ + • • ■ 

= n- P{n)A^^ 4 , 

M := SnS-^ = (e^^^") +(^i(n)A + ■ • • )n(e^(") +Mn)A + ■ • ■ 

= n + (^i(n)e-'^("+iU + -- - . 

Proposition [8] Let us take W of Theorem [5] and consider Qki ■= W^^CkiW which satisfy [8fc;,A] 
[&ki, n] = and hence Qki do not depends on A nor on n. Now, 

Ek'k'Qkl — Sk'k'dk'l, QklEk'k' — Sik'Qkl => Qki = ^klEkl, "dki £ C, 
Ek'k'Qkl = Sk'k^k'l, QklEk'k' — Sik'Qkl ©fci — "^klEkl, "Ski & C. 

Thus, 

Cki = WOkiW-' = 5VKoi9fciSfciVKo-i5-i = ^kiL'"-'' e^^^^'-'^"-'''^'^' {Eki + 2-) 

^dkt = l^ Cki = WEkiW-\ 



Acknowledgements 

The authors wish to thank the Spanish Ministerio de Ciencia e Innovacion, research projects FIS2005-00319 and 
FIS2008-00200, and acknowledge the support received from the European Science Foundation (ESF) and the 
activity entitled Methods of Integrable Systems, Geometry, Applied Mathematics (MISGAM). This paper was 
finished during the research visits of one of the authors (MM) to the Universite Catholique de Louvain and 
to the Scuola Internazionale Superiore di Studi ^?;anza<i/Intcrnational School for Advanced Studies (SISSA) in 
Trieste, MM wish to thanks Prof, van Moerbeke and Prof. Dubrovin for their warm hospitality, acknowledge 
economical support from MISGAM and SISSA and reckons different conversations with P. van Moerbeke, T. 
Grava, G. Carlet and M. Caffasso. 



25 



References 



[1] M. Adler and P. van Moerbeke, Comm. Math. Phys. 203 (1999) 185 . 
[2] M. Adler and P. van Moerbeke, Comm. Pure App. Math. 54 (2001) 153. 
[3] M. Adler, P. van Moerbeke, and P. Vanhaecke, Comm. Math. Phys. 286 (2009) 1. 
M.Adler. J. Delepine and P. van Moerbeke, Comm. Pure App. Math. 62 (2009) 334. 

E. Daems and A.B.J. Kuijlaars, J. of Approx. Theory 146 (2007) 91. 

[4] M. J. Bergvelt and A. P. E. ten Kroode, Pacific J. Math. 171 (1995) 23 . 
[5] A. Bottcher, M. Embree, and V. I. Sokolov, Lin. Alg. App. 343-344 (2002) 101. 
T. Strohmer, Lin. Alg. App. 343-344 (2002) 321. 

[6] M. Cafasso, Matrix biorthogonal polynomials on the unit circle and non-abelian Ablowitz-Ladik hierarchy. 
arXiv:0804.3572v2 [math.CA] 23 Apr 2008 

[7] G. Carlct, J. Phys. A: Math. Gen. 39 (2006) 9411. 

[8] E. Date, M. Jimbo, M. Kashiwara, and T. Miwa, J. Phys. Soc. Japan 40 (1981) 3806. 
[9] A. Doliwa and P.M. Santini, Phys. Lett. A 233 (1997) 365. 

M. Manas, A. Doliwa and P.M. Santini, Phys. Lett. A 232 (1997) 99. 

A. Doliwa, P.M. Santini and M. Manas, J. Math. Phys. 41 (2000) 944. 

[10] B. Eynard, An Introduction to Random Matrices, lectures given at Saclay, October 2000, http://www- 
spht.cea.fr/articles/t01/014/. 

[11] P. di Francesco, P. Ginsparg and Z. Zinn-Justin, Phys. Rept. 254 (1995) 1. 

[12] E. Getzler, The Toda conjecture, math. AG/0108108. Published in: Symplectic geometry and mirror sym- 
metry (Seoul, 2000), 5179, World Sci. Publishing, River Edge, NJ, 2001. 
G. Carlet, Theor. Math. Phys. 137 (2003) 1390. 

G. Carlet, B. Dubrovin, and Y. Zhang, Moscow Math. J. 4 (2004) 313. 

B. Dubrovin and Y. Zhang, Comm. Math. Phys. 250 (2004) 161. 
[13] J. Harnad. and A. Yu. Orlov, Theor. Math. Phys, 152 (2007) 1099 

[14] A. Gerasimov, A. Marshakov, A. Mironov, A. Morozov and A. Orlov, Nuc. Phys. B 357 (1991) 565. 
[15] V. G. Kac and J. W. van dc Lcur, J. Math. Phys. 44 (2003) 3245. 
[16] I. M. Krichcvcr, Comm. Pure. Appl. Math. 47 (1994) 437. 

[17] B. G. Konopelchenko and W. Oevel, "Matrix Sato Theory and Integrable Equations in 2+1 Dimensions" in 
Proceeding NEEDS'91, Baia Verde, Italy, June 1991. 

[18] E. J. Martinec, Comm. Math. Phys. 138 (1991) 437. 

[19] M. Mafias, L. Martinez Alonso, and E. Medina, J. Phys. A: Math. Gen. 35 (2002) 401. 
L. Martinez Alonso and M. Mafias, J. Math. Phys. 44 3294 (2003). 

F. Guil, M. Mafias, and L. Martfnez Alonso, J. Phys. A: Math. Gen. 36 (2003) 4047. 
F. Guil, M. Mafias, and L. Martfnez Alonso, J. Phys. A: Math. Gen. 36 (2003) 6457. 
M. Mafias, J. Phys. A: Math. Gen. 37 (2004) 9195. 

M. Mafias, J. Phys. A: Math. Gen. 37 (2004) 11191. 
[20] M. Mafias, L. Martfnez Alonso, and E. Medina, J. Phys.A: Math. Gen. 33 (2000) 2871. 

M. Mafias, L. Martfnez Alonso, and E. Medina, J. Phys.A: Math. Gen. 33 (2000) 7181. 
[21] L. Martfnez Alonso and E. Medina, J. Phys. A: Math. Gen. 40 (2007) 14223. 

[22] L. Martfnez Alonso and E. Medina, Multiple orthogonal polynomials, string equations and the large-n limit 
, arXiv: 0812. 3817. 

[23] A. Yu Orlov and E. I. Schulman, Lett. Math. Phys. 12 (1986) 171. 

[24] A. G. Reimann and M. A. Semenov-Tyan-Shanski, J. Math. Sci. 31 (1985) 3399. 

[25] K. Takasaki, Commun. Math. Phys. 181 (1996) 131. 

M. Mafias, E. Medina and L. Martfnez Alonso, J. Phys.A: Math. Gen. 39 (2006) 2349. 

L. Martfnez Alonso, E. Medina and M. Mafias, J. Math. Phys. 47 (2006) 083512. 



26 



[26] K. Takasaki and T. Takebe, Rev. Math. Phys. 7 (1995) 743. 

[27] K. Takasaki, "Dispersionless integrable hierarchies revisited", talk delivered at SISSA at September 2005 
(MISGAM program). 

[28] K. Takasaki and T. Takebe, Physica D 235 (2007) 109. 

[29] K. Takasaki and T. Takebe, Lowner equations, Hirota equations and reductions of universal Whitham 
hierachy arXiv:0808.1444 

[30] K. Ueno and K. Takasaki, Adv. Stud. Pure Math. 4 (1984) 1. 

[31] H. Widom, Trans. American Math. Soc. 121 (1966) 1. 

M. Van Barel, V. Ptak, Z. Vaverm, Lin. Alg. App. 332-334 (2001) 583. 
[32] P. B. Wiegmann and A. Zabrodin, Comm. Math. Phys. 213 (2000) 523. 

M. Minecv-Weinstcin, P.B. Wiegmann, and A. Zabrodin, Phys. Rev. Lett. 84 (2000) 5106 

I. M. Krichever, M. Mineev-Weinstein, P. B. Wiegmann, and A. Zabrodin, Physica D198 (2004) 1. 

L M. Krichever, A. Marshakov and A. Zabrodin, Comm. Math. Phys. 259 (2005) 1. 



27 



